The influence of structural short-range order on the phase diagrams 
of diluted FCC magnet with arbitrary spin 
and modified RKKY interaction 
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A diluted FCC magnet with modified long-range RKKY interaction and arbitrary Ising spin S is 
considered within two-sublattice model. In the molecular field approximation the Gibbs free-energy 
is derived, from which all magnetic thermodynamic properties can be self-consistently obtained. In 
particular, the phase diagrams are studied for different magnetic ion and free-charge concentration, 
the atomic short-range-order (Warren-Cowley) parameter being taken into account. 



I. INTRODUCTION 

The diluted magnetic systems constitute an important 
field of research in solid-state and statistical physics. The 
description of such alloys involves different challenges, 
starting from the question of the electronic structure of 
inhomogenous systems lacking translational symmetry 
through the calculation of exchange interactions between 
spins and finally to the problem of magnetic ordering. 
Another issue is the usual difference in the energy scales 
of magnetic and non-magnetic interactions between the 
ions in the alloy. Studies of such substances have recently 
attracted attention owing to the interest in diluted mag- 
netic semiconductors (DMS), important from the spin- 
tronics point of viewi 

The vital importance of disorder in description of DMS 
has recently been emphasised by numerous theoretical 
works and stimulated by the experimental data showing 
the great sensitivity of the critical temperature and mag- 
netic ordering in Ga^Mni-^As samples to their treat- 
ment even though the concentration of magnetic impu- 
rities remains unchanged. ? This suggests that the occu- 
pation of the lattice sites by magnetic ions may not be 
completely random. 

For the case of this celebrated DMS, Ga^Mni-^As, 
the effects of magnetic impurity clustering on the critical 
temperature of a model DMS have been investigated; for 
instance, by means of MC method in the work of Priour 
and Das Sarma£, who found only a weak influence of 
magnetic ion aggregating on the phase transition. By 
contrast, Bouzerar et alA predicted a noticeable increase 
in Curie temperature of the clustered impurities system 
for Ga^Mni.^As and Ga^Mn^jN. The presence of clus- 
tering in DMS has been supported theoretically, for ex- 
ample, by the first-principle calculations of various alloy 
characteristics performed by Drchal et al. 5 , Kudrnovsky 
et al£ or Raebiger et a/£ On the other hand, numerous 
existing works use different approaches to the clustering 
in order to show that it decreases the critical temperature 
in Ga-based DMS£&±&ii 

Various variants of ab initio methods, yielding the 
crystalline, electronic and magnetic structure, are be- 
lieved to provide reliable estimations of critical temper- 



atures for specific substances (for example see Ref. Qjl) . 
At the same time, however, it is difficult to extract from 
those methods the systematic analysis of the effect of 
particular factors on the final parameters of the system. 
Here we see the importance of simplified, schematic mod- 
els which could be extensively analysed, providing some 
detailed insight, for instance, into the general features 
of the influence of the disorder on the magnetic proper- 
ties. Our aim is to analyze a model diluted magnet with 
modified RKKY interaction, focusing on the importance 
of structural correlations in distribution of magnetic mo- 
ments on the lattice. We improve the virtual crystal ap- 
proximation by taking into consideration the correlations 
of pairs and study its effect on the magnetic ground state 
phase diagrams as well as the critical temperature. More- 
over, we introduce a measure of magnetic frustration and 
discuss its sensivity to structural correlations. 

The paper is organized as follows: In Ilnd section the 
theoretical model is described in detail by a statistical- 
thermodynamical method. In particular, the analytical 
expressions for the Gibbs energy and phase transition 
temperature are derived. On this basis the numerical cal- 
culations are carried out in Illrd section and the results 
are presented in figures. A recently developed approach 13 
of numerical summation over arbitrary large number of 
co-ordination zones is adopted. Discussion of the results 
in Illrd sec. is focused on the influence of atomic dilution 
and short-range-order on the magnetic phase diagrams. 
Finally, in IVth section some conclusions are drawn. The 
paper is supplemented by Appendix containing a method 
of configurational averaging in the pair approximation. 



II. THEORETICAL MODEL 

We will consider the Ising-type Hamiltonian on the di- 
luted FCC lattice containing antifcrromagnetic nearest 
neighbour (NN) interaction as well as the long-range in- 
direct interaction of the RKKY kind. In order to describe 
various antiferromagnetic structures a model of two in- 
terpenetrating sublattices (a, b) is adopted. The Hamil- 



2 



tonian can be written in the form: 
7~L = — E^ SiSj — E^ 



(1) 

where S" — —S, ... ,+S is the Ising spin of arbitrary 
magnitude S situated in i-th lattice site and belonging 
to the sublattice a (a = a,b). In Eq. (JTJ h — —g cS [iBH z 
corresponds to the external magnetic field H z oriented 
in z-direction, whereas g eS is the effective gyromagnetic 
factor, which for the case of RKKY interaction has 
been introduced in Ref. 0. The occupation operators 
& = (0, 1) describe the magnetic dilution. Namely, 
£i = corresponds to the magnetic vacancy in the i-th 
lattice site, & = 1 describes the state when i-th lattice 
site is occupied by the spin Sf. These operators are 
subject to configurational averaging (...) r , which we 
assume to be independent of the magnetic structure. 
On the other hand, the spin operators Sf, for a given 
atomic configuration, are subject to thermal averaging 
(...) only. 

The configurational averaging of single-site occupation 
operators can be conveniently described by introducing a 
parameter n = (£i) r , where n is a concentration of mag- 
netic component. The parameter n can be regarded as 
a quotient of the number of magnetic atoms to the total 
number of lattice sites. For simplicity, we will further 
assume that n is equal for both sublattices (a = a, b). In 
turn, the configurational averaging of -pairs leads 
to the expression: 



(2) 



where & = {£i) r + 8t;i, and A& = (8£i5£j) r is a fluctuation 
of the occupation numbers, which is characteristic for 
the fc-th coordination zone. As has been shown in the 
Appendix, the fluctuations must obey the sum rule: 



E Zk ^ k 







(3) 



where Zk is the total number of lattice sites on the fc- 
th co-ordination zone where the fluctuation takes 
place. These fluctuations are connected with the Warren- 
Cowley (W-C) short-range-order (SRO) parameter au by 
the relationship: 
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(4) 



A detailed analysis of the physical range of W-C pa- 
rameter resulting from the pair probability distribution 
P ^ s presented in the Appendix. 



As far as the thermal averaging of the spin operators 
is concerned, we will adopt the simplest molecular field 
approximation (MFA), with the decoupling relation 

(S?S^ « m a mP where a (or f3)= a,b and m a = (5f> 

denotes a-sublattice magnetization. The MFA is justified 
both by the presence of long-range interaction (for the 
infinite interaction range it becomes an exact method) 
and by the sublattice model of antifcrromagnetism 
according to the idea of Neel. As far as we know / 5 ' 16 ' 17 
in the FCC structure, apart from the ferromagnetic 
(F) and paramagnetic (P) phases, nothing impedes the 
notion that different antiferromagnetic orderings exist. 
The most known seem to be the antiferromagnetic 1st 
kind (AF1), antiferromagnetic 1st kind improved (AF1I) 
and 2nd kind antiferromagnetic (AF2) orderings. 



Within MFA the magnetic enthalpy can be found by 
the configurational and thermal averaging of the Hamil- 
tonian (JTJ. The result is: 



H=((H)) r = 



(m a f+ (m 6 ) 2 



N 



fe) TO TO 



— n (m a + m b ) h 



(5) 



By N we denote the total number of lattice sites, 
whereas the summation upon fc is performed over all 
co-ordination zones centered at the arbitrary lattice site. 
In Eq. © z\ ] ( zl l ) arc the number of lattice sites on 
the fc-th co-ordination zone, whose spins (if occupied) 
are oriented parallelly (antiparallelly) to the central 
spin. Thus, z k is the co-ordination number at the fc-th 
zone formed from lattice sites belonging to the same 
magnetic sublattice as the central spin, while zV' is the 
co-ordination number formed from lattice sites belonging 
to different sublattice. Those numbers satisfy the condi- 

T t T I 

tion z k + z k = Zk and their disribution upon fc depends 
on the type of magnetic ordering (F, AF1, AF1I or AF2). 

The exchange integral Jk in Eq. |5| for a given co- 
ordination zone fc is basically the RKKY long-range inter- 
action, with the exception of the first co-ordination zone, 
where we additionally include the antiferromagnetic su- 
perexchange interaction J AF < 0. This kind of inter- 



,18,19 



con- 



action has been introduced in several paper: 
corning diluted magnetic semiconductors (DMS). Thus, 
we assume that J\ — J AF + J 1 RKKY for fc = 1 and 
J k = JRKKY for k = 2,3,..., where the RKKY inter- 
action is given by the expression)^ 



J, 



RKKY 



C (k F a)' 



sin (2fcprfe) — 2fcprfe cos (2fcprfc) 
(2fc F r fc ) 4 



-T-fc/A 



(6) 



In Eq. © a is the lattice constant and rj, stands 
for the radius of the fc-th co-ordination zone. The 
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Fermi wavevector kp for FCC structure takes the form 

1 /3 

kp = (l27r 2 7i) / a, where we assume that each occupied 
lattice site yields one charge carrier to the conduction 
or valence band. Thus, we assume that the free carriers 
concentration is n, i.e., the same as the concentration of 
magnetic atoms in the FCC lattice. The energy constant 
C in Eq. ([6]) can be treated as the unit energy, both for 
the exchange integral and for the k B T scale as well. As 
far as the exponential factor in Eq. ((6J, containing the 
damping parameter A, is concerned, such a term has 
been introduced by Mattis 21 in order to account for the 
charge carriers localization in disordered systems. It 
follows from the literature ) 3 ' 19 ! 22 ' 23 that such localization 
takes place in some DMS, such as Gai-^Mn^As. The 
above description of J k for k = 1,2,... defines the 
so-called modified RKKY interaction, the model which 
has successfully been used in several papers concerning 
DMS systems^ 

The magnetic enthalpy (O allows the studies of the 
ground-state phase diagrams (for T — > 0) when a perfect 
spin alignment is assumed (for instance, m a — S and 
m b = — S for the antiferromagnetic phase). By compar- 
ison of the enthalpy values for different magnetic phases 
(F, P, AF1, AF1I and AF2) the stability areas for each 
phase can be established from the minimum condition. 
The magnetic phase diagram can be obtained vs. concen- 
tration n, fluctuation distribution A k and external field 
h, for given parameters J AF and A characterizing the 
modified RKKY interaction. For instance, for the an- 
tiferromagnetic phases in the ground state the enthalpy 
per lattice site is given by the expression: 



H 

N 



k 



(l + a k )(zl r -zl L )j k (7) 



where a k is W-C parameter for the fc-th co-ordination 
zone and zt {zl ) depend on the type of the antiferro- 
magnetic phase. On the other hand, for the ferromag- 
netic phase in the ground state (when m a — m b = S) 
the corresponding formula for the enthalpy reads: 



H 

N 



i 7J 2 S 2 y^(l + a k ) z k J k + - nSh 

k 



(8) 



In turn, for the paramagnetic phase (with m a = 0) 
we assume H/N = 0. The numerical calculations of 
the ground state phase diagrams with the structural 
ordering taken into account will be presented in the next 
section. 

As far as the temperature studies are concerned, one 
has to consider not only the temperature dependencies 
of the sublattice magnetizations m a (a — a, b) but also 
the magnetic entropy. In the MFA method a unified ap- 
proach can be suggested, based on the single-site density 
matrix 



cxp \/3 (A a + h) S[ 



(9) 



where j3 — 1/ksT, and A" is a variational parameter of 
the molecular field acting on a-sublattice. The single-site 
partition function for a-sublattice is defined by the 
formula: 

s 

Z? = Tr 4 {exp \p (A Q + h) Sf]} = ^ exp [/3 (A Q + h) I] 

i=-s 

(10) 

Then, the total partition function Z in the MFA is given 
by the product: Z = (Z?Z*l) Nn/2 , where Nn/2 is the 
number of lattice sites occupied within one sublattice. 
With the use of single-site density matrix @ the various 
thermal mean values can be calculated. For instance, the 
magnetization of the occupied site on a-sublattice is: 

m a = Tr, [Sfpf] = i-Tr^Sf exp [3 (A<* + h) S? }} 



which leads to the general formula 

m a = SB(S/3(A a + h)) 



(11) 



(12) 



where a = a, b and SB (x) is the Brillouin function for 
an arbitrary spin S: 



^(S.) = ^±W^±i. 



icoth(|) (13) 



Due to the factorization in MFA the total entropy a 
can be presented as a sum of single-site entropies for both 
sublattices: 



N 



® = "o" 1^ CT 1 



(14) 



a— a, b 



The single-site entropy a" for the occupied site on a- 
sublattice is given by the thermal mean value 

< = -k B (Inpf) = -fc B T ri (pflnpf) (15) 

Making use of Eq. (J9j> , the single-site entropy af can be 
presented in the form: 

af=~ (A" + h)m a + k B lnZ? (16) 

hence the total entropy (fTl| can be presented as: 
N 1 



a = — n — 
2 T 



(A a + h) m a - (A b + h) m b 



k B TlnZ? + fc B TlnZ 1 h ] 



(17) 



Having calculated the entropy (fT7|) and the enthalpy 
(JSJ) , the Gibbs free-energy can be found in MFA from the 
thermodynamic formula: 



G = H -oT 



(18) 
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However, with a view to constructing the phase diagrams, 
we are interested in the Gibbs energy per lattice site, i.e., 
in the chemical potential fj,, which is given by: 



G 
N 



1 



$>z£ T (l + a fc ) [(m a ) 2 + ( m b Y 

k 

\n 2 Y,Jkzl l {l + a k ) m a m b 



+| (A a m a + A b m b ) - ^k B T (lnZ? + lnZ?) 

(19) 



The molecular field (variational) parameters A Q appear- 
ing in Eqs. (fT9|) . (fT2|) and (fTOf can be determined from 
the necessary extremum conditions: 







(20) 



(for a = a,b) which lead to the following expressions: 



A a = ]Tj fc (l + a fe ) ( 



z[}m a + z\ i m h 



and 



A 6 = ^J fc (l + a fc ) (4 T m b + 



z\ l m a 



(21) 



(22) 



Now, with the help of Eqs. (|2Tj) and (|22|) the chemical 
potential (|19p for the thermodynamical equilibrium is ob- 
tained in the final form: 



1 



; ^J fe 4 T (l + a fc ) \(m a f + (m b ) 



(23) 



together with Zf given by Eq. (fTU|) and m Q as a solution 
ofEq. (HU). 



From this point on, the expression (f2"3"|) for the chem- 
ical potential allows the self-consistent studies of all the 
thermodynamic properties in MFA. Let us notice first 
that the necessary equilibrium conditions for the chem- 
ical potential (f2T)]) are also satisfied with respect to the 
sublattice magnetization: dfj,/dm a = (a = a, 6). More- 
over, the mean magnetization per one lattice site m can 
be derived alternatively to Eqs. (fTTl [T2")) merely by dif- 
ferentiation of the chemical potential over the external 
field: 



dfi\ 1 , a 
— — = —n to 
dh) T 2 v 



(24) 



Analogously, by differentiation of the chemical potential 
over temperature the mean entropy per lattice site can 



be calculated, yielding the same form as Eqs. (|14I17|) : 



N 



df 



\n {<% + <#) 



(25) 



Consequently, other thermodynamic properties such as 
the magnetic susceptibility, or magnetic contribution to 
the specific heat can be calculated as the second-order 
derivatives of the chemical potential (j2"3"]) . 

The critical temperature of the second-order (contin- 
uous) phase transitions can be obtained from the lin- 
earization of Eq. (fT2"j) for h — and m a — » 0. Mak- 
ing use of the linear expansion for the Brillouin function 

SB (Sx) S(S+l) x/3 we obtain from Eq. (JT2J: 



S{S + l) 



(26) 



(a = a, 6), where /3 C = l/fceTc and T c is the critical tem- 
perature. Now substituting A a from Eqs. (|2ip and l]22p 
into (|26p . we obtain a set of two linear, homogeneous 
equations for m a — * in the vicinity of T c . Next, by 
setting the determinant to be equal zero, the phase tran- 
sition temperature T c is derived in the following form: 



S(S + 1) 



,^J fc (1 + a k ) 



z k 31 z k 



(27) 



Eq. ((27|) is a generalization of MFA result for the 
long-range RKKY interaction with the structural 
clustering taken into account. The solution with "+" 
corresponds to the Curie temperature and is applicable 
to ferromagnetic phase transition, whereas the solution 
with "— " corresponds to the Neel temperature for 
the antiferromagnetic (AF1, AF1I, and AF2) phase 
boundaries. The temperature phase diagrams based on 
Eq. (|2"7| will be calculated in the next section. 



III. NUMERICAL RESULTS AND DISCUSSION 

The numerical studies have been carried out for the 
model FCC structure with dilution and the modified 
RKKY interaction taken into account. The external 
magnetic field h = was set. According to the theo- 
retical considerations, when the dilution is not random 
it cannot be desribed by only a single variable 77, but 
should be characterized by the set of Warren-Cowley 
SRO parameters. These parameters, a k , which are 
incorporated into the analytical formulas in previous 
section, fulfil the sum rule (Eq. [3]) and are treated as 
independent of the magnetic structure. However, for 
the simplicity of numerical calculations, we will further 
assume that only ot\ and a-i parameters are different 
from zero, while a k = for k > 2. Such assumption 
reflects the empirical fact that the structural correlations 
for two first coordination zones are the most important 
factors. In this way only a\ becomes an independent 
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FIG. 1: Ground-state magnetic phase diagrams in the (n,Qi)-plane for different values of parameters modifying the RKKY 
interaction: A — > oo, J AF = (a), A -> oo, J AF /C = -0.5 (b), A/a = f .0, J AF = (c), A/a = 1.0, J AF /C = -0.5 (d). By U 
we denote an unphysical area for ai-parameter. 



SRO parameter of the theory (with some constraints 
imposed, as discussed in the Appendix (see IA. 161 and 
Fig. 6)), whereas oti is determined from the sum rule: 
a 2 = - (zi/z 2 ) ai. 

All numerical calculations involving summation 
over co-ordination zones k have been performed up 
to fc max = 41253, which corresponds to the radius 
of interaction rfc max /a = 150 in the FCC structure. 
The co-ordination numbers zjj (zt ) for a particular 
magnetic phase have been computed by a computer 
program analysing the lattice. For such large range 
of interaction a perfect numerical convergency of the 
summation has been achieved. Hence, for the results 
presented below, the numerical errors are negligible. In 
particular, in case of the ground-state phase diagrams, 
where no analytical approximations have been made, 
the presented results can be regarded as near-exact. 



A. Ground-state phase diagrams 

The ground-state phase diagrams have been computed 
on the basis of Eqs. and © for h = 0. The stability 



regions for the phases (F, P, AF1, AF1I and AF2), where 
the enthalpy of each phase is in minimum, have been 
found in (n, ai)-plane, and the results are presented 
in Fig. [1] The (a)-(d) parts of Fig. [1] correspond to 
the different parameters of modified RKKY interaction, 
as indicated in the figure caption. In all parts (a)-(d) 
the unphysical (U) area of Warren-Cowley parameter 
oil vs. charge carriers concentration n is depicted and 
delimitated by the dashed lines. The origin of that area 
was discussed in the Appendix. The horizontal dashed 
line for a>i — corresponds to te absence of SRO. The 
ocurrence of phases and their sequence for ot\ = upon 
n is in agreement with those studied in Refs. Il3ll24l . 
where the SRO was not taken into account. 

Fig. 1(a) is prepared for A — > oo and J AF = 0, i.e., for 
the "pure" RKKY interaction. The ferromagnetic (F) 
phase is present for the low concentration n only, but for 
all values of a\. The positive a\ values (with structural 
clustering) enlarges the stability of the F phase. On the 
other hand, for a± < the antiferromagnetic phases 
become more favourable, and even the paramagnetic (P) 
phase can occur in some small restricted area. There 
are three stability regions for AFl-phase in Fig. 1(a). 
The rest of phases occur in the single regions only. In 
Fig. 1(b) the antiferromagnetic NN interaction, with 
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FIG. 2: MFA reduced critical temperature vs. charge carriers concentration, for different values of parameters modifying the 
RKKY interaction: A — > oo, J AF = (a), A -> oo, J AF /C = -0.5 (b), A/a = 1.0, J AF = (c), A/a = 1.0, J AF /C = -0.5 (d) 
and for different structural correlations: ot\ — (solid line), Qi = —1.0 (dashed line) and «i = 0.5 (dashed-dotted line). 



the value J AF /C = —0.5, is included. As a result the 
F-phase area is strongly reduced, whereas the antiferro- 
magnetic phases are extended. There are four stability 
regions for AFl-phase, and two for each AF1I and AF2. 
The new areas (AF1, AF1I and AF2) appear for the 
lowest values of n. Moreover, the P-phase becomes 
stable over all values of a\ in some narrow range of 
concentrations n. In Fig. 1(c) we assume J AF = and 
A/a = 1, which represents a relatively strong charge 
carriers localization. The topology of this diagram is 
somewhat similar to Fig. 1(a), with some differences 
in the shapes of the particular lines. Two areas of 
antifcrromagnetic AFl-phases (one for strong negative 
and other for strong positive ai-parameters) are smaller 
than the corresponding areas in Fig. 1(a), i.e., without 
the RKKY damping. In the last part (Fig. 1(d)) both 
non-zero values of J AF /C = —0.5 and A/a = 1 were 
assumed. As a result, the F and P phases have been 
reduced to the very small areas. It is seen that for 
the positive structural clustering, with oil > 0, only 
the antifcrromagnetic phases are in favour in Fig. 1(d). 
There are two AF1 and also two AF1I areas which 
are stable for these phases in that figure. The AF1 
region has remarkably increased, whereas the AF2 phase 
has been limited mainly to the negative SRO parameters. 



From the analysis of Figs. l(a)-(d) some general 
conclusions can be drawn. The SRO parameter exerts 
a strong influence on the phase diagram, but not for 
all values of n on the equal footing. The greatest 
changes are observed for small n values, i.e., for strong 
dilution. Both J AF and A parameters, which modify 
the RKKY interaction, have also important meaning, 
which is manifested in Fig. 1 by remarkable divergences 
between various parts ((a)-(d)) of the figure. All the 
phase boundaries presented in Fig.l constitute the 
discontinuous (1st order) phase transitions. 



B. Critical temperature 

The transition temperature from spontaneously or- 
dered to disordered state is calculated on the basis of 
Eq. (|27|) . In connection with the ground-state phase 
diagrams, the most important changes are predicted 
for small n- values; therefore we will present the results 
of critical temperature calculations in the range of 
n < 0.25. In Fig. 2, in the parts (a)-(d), the reduced 
critical temperature T c / [nS(S + 1)C] vs. n is presented, 
for the same parameters of modified RKKY interaction 
as in Fig. 1 in the parts (a)-(d), respectively. The 
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reduced temperature provides the results indepen- 
dent of the spin magnitude S and shows a non-linear 
dependency of T c upon n. On each phase diagram 
(2(a)-2(d)) the three curves are presented: solid, dashed 
and dashed-dotted, for ot\ = 0, —1 and 0.5, respectively. 
Thus, the temperature phase diagrams in Fig. 2 are 
prepared as the horizontal cross-sections for the cor- 
responding ground-state phase diagrams presented in 
Fig. 1. The values of cti have been chosen on the basis 
that enables to compare the calculations for the extremal 
SRO-parameters (c*i = — 1 and a\ = 0.5) with that for 
the absence of SRO (ct\ = 0). The critical temperatures 
for each specific phase in Fig. 2 are pointed by the 
arrows with the phase symbols. The vertical (dotted) 
lines indicate the discontinuous phase boundaries, from 
P-phase down to the ground state, and their positions 
upon n are in agreement with Fig. 1. 

Some interesting features observable in Fig. 2 should 
be emphasized. In particular, in Fig. 2(a), for the un- 
modified RKKY interaction, it is seen that the positive 
clustering increases the Curie temperature and makes the 
region of F-phase wider. On the other hand, the negative 
structural correlations reduce both the Curie tempera- 
tures and the width of the ferromagnetic region. For the 
negative SRO parameters, with increase of n a relatively 
high Neel temperatures for AF2 phase can be achieved. 

For the antiferromagnetic NN interaction with the 
value J AF /C = —0.5 (Fig. 2(b)) the situation seems 
somewhat more complicated. The ferromagnetic phase 
for a± = and a\ = 0.5 is weaker than for the corre- 
sponding cases in Fig. 2(a); however, for a,\ = —1 it 
has a very similar phase boundary. The characteristic 
features are three paramagnetic gaps (one for each 
curve) occurring in accordance with Fig. 1(b). Also, 
for the small values of n the new phases AF1 and 
AF1I are perceptible, with their Neel temperatures 
decreasing with increase of magnetic impurity (and 
charge) concentration n. In Fig. 2(c) the influence of 
the charge carriers localization (with A/a = 1) on the 
Curie temperature is demonstrated. When compared 
with Fig. 2(a), a remarkable reduction of all critical 
temperatures is apparent although the influence of 
SRO on the phase boundary shapes is comparable for 
both figures (1(a) and 2(c)) over the presented range 
of n. In Fig. 2(d), when both non-zero parameters 
(J AF /C = —0.5 and A/a = 1) are taken into account, 
the ferromagnetic phase is present only for a% = —1. 
For all curves in Fig. 2(d) a minimum of T c is present in 
the region of concentrations for n « 0.07 — 0.15. 




FIG. 3: MFA critical temperature in the vicinity of a se- 
lected triple point, (n,ai) = (0.077,-0.25) from Fig. 1(b), 
vs. Warren-Cowley parameter ai, for various charge carriers 
concentrations n. 



0.98 



- 0.96 



0.94 




0.92 



FIG. 4: Degree of frustration / in the vicinity of a selected 
triple point with the same parameters as in Fig. 3, vs. Warren- 
Cowley parameter ai, for various charge carriers concentra- 
tions n. 



transition lines merge with another 1st order phase 
transition line are usually called the bicritical points. It 
can be noted that with the change of concentration n 
the phase transition temperatures can either increase or 
decrease, and the changes are strongly non-linear. This 
effect depends not only on the range of concentration n 
but on the SRO parameter as well. 



All the phase transitions from ordered to paramagnetic 
state, shown in Fig. 2, are continuous (2nd order) ones. 
However, for the vertical (dotted) lines, the perpendicu- 
lar phase transitions (for constant temperature) between 
different ordered states are of discontinuous character, 
starting from P-phase down to the ground state. The 
points of junction, where two different 2nd order phase 



C. Magnetic frustration 

Magnetic frustration is unavoidable for systems with 
the long-range oscillatory interaction. The presence of 
frustration leads to the increase of magnetic energy and 
decrease of the critical temperature. In order to define 
quantitatively the degree of frustration, we introduce an 
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FIG. 5: Distribution of /-parameter over (n, ai)-plane for J F /C = and A — > oo. The contour lines connect points with the 
same value / of magnetic frustration. Note the differences in / scale in the diagrams. 



'ideal" ground-state energy Eq in the form 

N 



E = -J2\ J ( r a)\ 

(i,3) 



^z k \J(r k )\. (28) 



Eq is the lowest magnetic energy which would be achieved 
in the hypothetical case, when no frustrations were 
present in the system. However, the real internal energy 
in the frustrated ground state, E g = (H) h _ , is given 
by the formulas (J7J) and ([5]) for the external field h = 0. 
Thus, we can define the degree of frustration / in the 
ground state as follows: 



/ 



E„ — Eq 



E Q 



(29) 



The value of / comes from the range (0; 1), where 
/ = is for the "ideal" unfrustrated system, and / = 1 
stands for the paramagnetic phase with E g = 0. 

In order to demonstrate a correlation of the degree of 
frustration with the phase transition temperature, first 
in Fig. 3 we present the four reduced T c -curves for dif- 
ferent n parameters, when J AF jC = —0.5 and A — > oo. 
The choice of n-parameters (n = 0.06, 0.077, 0.09 
and 0.11) enables the vertical scan of Fig. 1(b) upon 
ot\, in the vicinity of the triple point which has the 
coordinates (n, a\) — (0.077,-0.25). In this triple 
point in Fig. 1(b) the three phases: F, AF1 and P do 
coexist. The temperature studies in Fig. 3 show both 
Neel and Curie temperatures, as well as the existence 



of paramagnetic (P) phase, where T c = 0. On the 
other hand, in Fig. 4 the degree of magnetic frustration 
/ in the ground state is presented, for the same set 
of parameters (J AF , A, n), and abscissa axis ai, as in 
Fig. 3. By comparison of Figs. 3 and 4 it is seen that 
the increase (decrease) of the critical temperature is 
accompanied by the decrease (increase) of the degree 
of frustration in the ground state, respectively. Both 
T c and / are the linear functions of ot\, showing a 
non-smooth behaviour upon ai at the 1st order phase 
transitions. A relatively high degree of frustration in 
Fig. 4 is remarkable, with its maximum in the vicinity 
of the triple point. It is worth noticing that at the triple 
point the chemical potentials of all coexisting phases 
are equal, and, because of the presence of P-phase, / = 1. 

In order to see how magnetic frustrations are dis- 
tributed over (n, ai)-plane, in Fig. 5 the /-parameter 
is presented in the form of iso-/ contour lines. The 
interaction parameters in Fig. 5 amount to (J AF = 0, 
A -> oo) (a), (J AF /C = -0.5, A -> oo) (b), (J AF = 0, 
A/a = 1) (c) and (J AF /C = -0.5, A/a = 1) (d), i.e., 
are the same as in Figs. l(a)-(d), respectively. The 
unphysical range of SRO-parameter is again denoted by 
U. The iso-/ contour lines connect the points with the 
same value of the frustration parameter /, according 
to the figures legend. For two figures ((a) and (b)) 
a relatively high values of / can be noticed. On the 
other hand, in Fig. 5(c) a remarkable reduction of 
the /-parameter can be observed, when compared 
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with Fig. 5(a). It is worthy to remind that respective 
phase diagrams (Fig. 1(a) and Fig. 1(c)) have a similar 
topology and have been prepared with the difference in 
A-parameter only. Thus, it is hard to escape the obvious 
conclusion that the damping of RKKY interaction 
reduces the degree of frustration. In contrast, a presence 
of NN antiferromagnetic interaction on the FCC lattice 
increases the frustrations, which are then unavoidable 
even for the neighbouring spins. A small white areas 
in Figs. 5(a)-(d) with / = 1 correspond to the para- 
magnetic phases as seen in Figs. l(a)-(d), respectively. 
However, the remaining contour lines of / do not reflect 
the shapes of the ground-state phase boundaries. The 
distribution of /-parameter vs. a\ shows that the SRO 
has a strong influence on the frustration. However, this 
effect depends simultaneously on the magnetic dilution n. 



IV. CONCLUSION 



The magnetic phase diagrams for the diluted FCC 
lattice with modified RKKY interaction have been 
studied, with SRO parameter taken into account. The 
formalism necessary for structural averaging in the pair 
approximation has been presented in the Appendix. As 
a result, the physical range of Warren-Cowley parameter 
has been established. In the theoretical part, the general 
statistical-thermodynamical method has been presented. 
On the basis of the above formalism, the phase diagrams 
and other magnetic properties can be studied in the 
external field h. The method can also be adopted 
for other kinds of localized spins Hamiltonians. The 
numerical calculations show that the SRO parameter 
ai has a remarkable influence on the phase diagrams, 
both in the ground state and for the temperature depen- 
dency, especially for low n. Both positive and negative 
structural ordering parameters have been considered. 
A comparison of the results with those for the virtual 
crystal approximation can be easily made by assuming 
on = 0. The temperature phase diagrams contain the 
vertical phase boundaries, which represent 1st order 
phase transitions for T > 0. The diluted system with 
RKKY interaction is mostly a very frustrated one, and 
the frustration parameter / depends both on the SRO 
and the parameters of modified interaction. The inverse 
correlation between the degree of frustration and the 
critical temperature has been established. One of the 
novelty of the presented method lies in incorporating 
the Warren-Cowley parameter into the magnetic model, 
which can be easily extended to other systems, where 
the structural clustering can play a role in magnetic 
phenomena. 
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APPENDIX: ATOMIC SHORT RANGE ORDER 
IN DILUTED SYSTEMS 

We consider a system of localized spins on the crys- 
talline lattice. In such a system, the spin dilution can 
be conveniently described by means of site occupation 
operators £^ for each site, possessing the eigenvalues 
(the corresponding i-th site is empty) and 1 (the i-th site 
is occupied by the magnetic impurity ion). The prob- 
ability distribution of the eigenvalues of the single-site 
occupation operator is as follows (see ex. Ref. \2w : 

p(&)=n*(&-l) + (l-n)«(£i). (A.l) 

where 8 (£, — 1) and <5(£i) are the Kronecker's deltas. It 
can be verified that =o \ P (£«) = 1> thus the distribu- 
tion is nomalized. 

The configurational average of the operator equals 

so that it is the occupation probability for a single 
site (i.e. average number of magnetic impurities per 
site). The average (£j) r is independent on the i-th site 
location, which comes from the fact that the distribution 
(jA.ljl is valid for the whole sample. 

Let us consider a general probability distribution of 
the occupation of a pair of sites i and j: 

P (6, &) = P%6 m S &) + p%6 &) 8 & - 1) 

+pp (£ - 1) 8 & ) + pl}8 (& - 1) 8 - 1) . 

(A.3) 

The number pf^ B is the probability of the event AB 
(A, B — 0,1) for the pair of lattice sites i and j. 
The normalization requires: 

t%+t%+P$+P& = l- ( A -4) 

The simplest approach to the description of the diluted 
system, so called Virtual Crystal Approximation (VCA) , 
relies on the assumption that the impurity ions are dis- 
tributed randomly in the lattice sites, thus they are mi- 
correlated. Then the occupations of the specified sites 
are statistically independent events and the probabil- 
ity distribution (|A.3[) takes a product form = 
P (6) P (&)> so = (1 - n) 2 , plj = p}° = n (1 - n) and 

1 1 2 

Ptj = nr. 
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In such a situation, the only parameter describing the 
distribution of magnetic impurities is n. 

However, VCA does not include the possible Short 
Range Ordering (SRO) in the diluted system, which takes 
place when the occupations of specific sites are not inde- 
pendent events, implying then p 7^ P (&) P ■ 
In particular when p(£,i,£,j) > P P we deal with 
clustering. The existence of SRO stems from the interac- 
tions between the impurity ions (especially the coulombic 
ones). 

The general form of the two-site probability distribu- 
tion (| A.3|) must be reducible to one-site distributions: 



E p(6 

6=0,1 



&)=p&) 



E p(&.&)=pte 
6=0,1 



(A.5a) 



(A.5b) 



The above conditions, together with Eq. IA.11 yield: 

„00 



P 



p^+pI 1 



'13 



pI'+pI 1 



y 1 «j 
Using Eq. IA.4I we obtain: 



= n 

= 1 -n 

= n 

= 1 - n. 



P 



00 



l-2n 



P 



11 



P 



01 



P 



10 



P 



11 



(A.6a) 



(A.6b) 



(A.7) 



Thus, the general two-site probability distribution can be 
parametrized independently by n and p}j . 

The average of the pair occupation operator reads: 



(HiWr = E ^iPte-^: 
6,6=0,1 



: Pi 



(A. 



thus it equals the probability of occupation of both sites 
i and j simultaneously. For VCA, the average factorizes: 

{CiQr = (&)r (Cj')r due t0 the P roduct form «f JJ (&, £,-)• 

Let us write the site occupation operator as: 



& = (&> r + <5&, 



(A.9) 



where (<5£,} r = to fulfill Eq |A.2l Thus, in a general case 
we obtain: 



(A.10) 



Let us denote (S£,iS£,j) r — A k . The above expression 
describes isotropic correlations, i.e. the value (£i£j) r de- 
pends solely on the distance r& between sites i and j (the 
j-th site lies on the &;-th coordination zone of the i-th 
site). In particular, we have Ak = in absence of SRO, 
for VCA. In the above notation the probabilities are: 



Pij 



= 1 - 2rH 
= n 2 + A k 



A k 



= n — n 2 — A k 
(A.ll) 



One can also write the conditional probabilities: 



Oil 

Pij 



vM 1 



1 - n 2 - A k 

1 - n 
i 2 + A, 



010 , n-n 
Pi = 1 - 



1-71 



A fc 



n 



(A.12) 



A I B 

where p-' is the conditional probability of an event 
B = 0, 1 for the site j under condition that an event 
A = 0, 1 occurred for the site i. It can be verified that 
the above probabilities obey Bayes theorem as well as 
complete probability theorem. 

The possible physical range of Afc follows from the 
interpretation of p?", p®j, p\® and p\j as probabilities, 
obeying < pfj B < 1. Using Eq. IA.11I we arrive at the 
set of inequalities: 



n + 2n 

^2 1 „ 



- 1 < A fc < ~n l + 2n 

- 1 < A fe < -n 2 + n 
-n 2 < A k < -n 2 + 1, 



(A.13) 



which must hold for every k. 

There exists another condition to impose on the pa- 
rameters Afc. Let us consider the operator of the total 
number of occupied pairs, 
age: 



2^,.,,, -i*..- with the aver - 



In the equation above, Zfc denotes the number of lattice 
sites on the fc-th co-ordination zone. If we deal with a 
constant number of impurity ions in the lattice, then the 
total number of occupied pairs is also constant and does 
not depend on the SRO existence. As mentioned above, 
Afc = in absence of SRO. Therefore, by comparison of 
Eq. IA. 141 taken for the absence of SRO and for arbitrary 
choice of impurity correlations we obtain a constraint: 



^z k A k 



0. 



(A.15) 



Let us restrict to the situation, when SRO is limited to 
nearest- and next-nearest neighbours of each site, i.e. 
when Ai ^ and A 2 ^ 0, whi le A fe = for k > 2. 
Notice that the constraint (|A.15l) excludes the possibil- 
ity that only a single value of Afc, for instance for n.n, 
is non-zero. Therefore, the situation we selected involves 
minimal number of Afc parameters. Since the constraint 
gives A2 = — (zi/ ' z%) A\ the distribution of impurities is 
then described by two independent numbers n and Ai. 

Writing the inequalities (|A.13j) for Ai and for A2 ex- 
pressed by Ai, we arrive at the allowed range of the pa- 
rameter Ai, which has to fulfill the following 12 inequal- 
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FIG. 6: The allowed range of Warren-Cowley parameter for 
the first co-ordination zone of the fee lattice. The dashed 
lines (12 in total) are the limits obtained from the inequalities 
(|A,16[) . The thick solid line encloses the allowed range of Qi 
for various n. The dotted line for ai = corresponds to 
vanishing SRO. 



ities: 



Z 2 I 2 

— (n 

z\ 

Z2 l 2 

— (n 

z\ 

z 2 / ; 

— (n 

z\ 



- 1 


< 




< 


-n 2 + 


2n 


- 1 


< 


Ai 


< 


-n 2 + 


n 


-n 2 


< 


Ai 


< 


-n 2 + 


1 


2n) 


< 


Ax 


< 


Z 2 ( 2 
Z\ 


-2n + l 


-n) 


< 


Ai 


< 


Z 2 I 2 
Z\ 


-n+l) 


-1) 


< 


Ai 


< 


z 2 2 

— n . 

Z\ 





(A.16) 



The existence of SRO in diluted systems can also be 
described conveniently by means of the Warren-Cowley 
(W-C) parameters aij,— 7 defined as: 



(^■) r -(£U£i), 



(A.17) 



These parameters can be equivalently given in the form 

Oiij = — ^1 — plj 1 /nj , where pJJ 1 = pjj jn is a conditional 

probability of occupying the j-th site if the z-th site is 
occupied. 

In our notation the W-C parameter for fc-th coordina- 
tion zone is given by: 



a k = A k /n 2 



(A.18) 



and the configurational average of pair occupation oper- 
ator reads: 



(tej) r = n 2 (l + a k ). 



(A.19) 



The allowed range of W-C parameter a± values for Ai 
fulfilling the inequalities (|A.16p for fee lattice is presented 
in the Fig. [SJ where it is bounded by thick solid lines. In 
this range the physically possible correlations are given 
by Eq. (|A.19[) . For other lattices characterized by specific 
sets of z k numbers the range of W-C parameter requires 
separate calculations based on the inequalities (|A.16[) . 
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